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Abstract 


Soil erosion is one of the most important environmental problems in the world, which has created 
many problems for human societies. Qualitative and quantitative assessment of erosion using 
experimental models for estimating erosion and sediment is one of the solutions through which soil 
erosion can be controlled to some extent and its amount can be minimized. The present study was 
conducted to compare EPM and RUSLE models in estimating erosion and sediment in 
Chamgardalan watershed of Ilam province. In this regard, the factors affecting erosion include rain 
erosivity, soil erodibility, topography, vegetation and soil conservation operations as well as 
observable erosion coefficients, rock and soil susceptibility to erosion, land use and average slope, 
used as a percentage used in the mentioned models and finally the amount of erosion and sediment 
produced was calculated according to the proposed models and compared with the amount observed 
in the sedimentation station of Chamgardalan watershed. The results show that the estimated 
sediment rate of RUSLE model was 16.37 tons per hectare per year and in the EPM method was 
21.49 tons per hectare per year, which is closer to the estimated value of RUSLE. Observed sediment 
of 16.58 tons per hectare per year, the efficiency of RUSLE model in estimating sediment in the 
study area is higher. Finally, it is recommended to use the EPM model for the initial phases of 
justification, overview and identification studies and RUSLE method for detailed and more accurate 
studies of erosion and sedimentation of the watershed. 
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حکیده 


ارزیابی کیفی و کمی فرسایش با استفاده از مدل‌های تجربی برآورد فرسایش و رسوب. یکی از راهکارهایی است که از طریق آن 
می‌توان فرسایش خاک را تا حدودی کنترل و مقدار آن را به حداقل رساند. پژوهش حاضر به‌منظور مقایسه دو مدل EPM‏ و 
RUSLE‏ در برآورد فرسایش و رسوب در حوضه چم‌گردلان استان ایلام انجام شده است. در این راستاء از عوامل مؤثر بر 
فرسایش شامل فرسایندگی Obl‏ فرسایش‌پذیری SE‏ توپوگرافی. پوشش گیاهی و عملیات حفاظت خاک و همچنین ضرایب 
فرسایش مشاهده‌ای» حساسیت سنگ و خاک به فرسایش» استفاده از زمين پا کاربری زمین» شیب متوسط زمین به درصد مورد 
استفاده در مدل‌های یاد شده استفاده شده و درنهایت مقدار فرسایش و رسوب تولید شده با توجه به مدل‌های ارائه شده محاسبه 
و با مقدار مشاهده شده در ایستگاه رسوب‌سنجی حوضه چم‌گردلان مورد مقایسه قرار گرفت. نتایج به دست آمده نشان می‌دهد 
oli‏ رسوب برآوردی مدل RUSLE‏ ۱۱/۳۷ تن در هکتار در سال و در روش YVES EPM‏ تن در هکتار در سال برآورد 
گردید که با توجه به نزدیک‌تر بودن مقدار برآوردی RUSLE‏ به مقدار رسوب مشاهده‌ای ۱۳/۵۸ تن در هکتار در Sle‏ کارایی 
مدل RUSLE‏ در برآورد رسوب منطقه موردمطالعه بیش‌تر است. درنهایت توصیه می‌شود. برای فازهای Syl‏ مطالعات 
توجیهی. اجمالی و شناسایی از مدل EPM‏ و برای مطالعات تفصیلی و دقیق‌تر برآورد فرسایش و رسوب حوضه از روش 
RUSLE‏ بهره گرفت. 


کلیدواژه‌ها: فرسایش خاک. برآورد رسوب. مدل RUSLE‏ » مدل .EPM‏ حوضه چم گردلان. استان ایلام 
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معد مه 

خاک یکی از مهم‌ترین عناصر طبیعت و با ارزش‌ترین بستر تولید برای بشر است. مهم‌ترین عاملی که وجود 
خاک را تهدید می‌کند. فرسایش است. امروزه با پیشرفت دانش و فن و با ماشینی شدن جهان روبرو هستیم 
که بیش از هر زمان دیگر منابع آب‌وخاک را تهدید می‌کنند. به‌صورتی که هرساله شاهد تولید میلیون‌ها تن 
رسوب در سراسر دنیا هستیم )1995 (Rafahi & Nemeti,‏ که فرسایش خاک بر اثر عامل آب بیش از ۵7 
درصد آن را شامل می‌شود )2001 (Elirehema,‏ روش‌های اندازه‌گیری شدت فرسایش در حوضه‌های آبخیز 
شامل روش‌های تجربی» روش‌های ریاضی و نمونه‌برداری رسوب است „Galal, Hegazy, & Alghailani,‏ 
(Bahrawi, Elhag, Aldhebiani 2016)‏ روش‌های تجربی برآورد فرسایش و رسوب در bul‏ شامل معادله 
جهانی تلفات خاک (Wischmeier & Smith, 1978 ) (USLE)‏ معادله اصلاحی تلفات خاک ‘(MUSLE)‏ 
L (Ferro, Giordano, & Lovino, 1991)‏ معادله تغییر cub act,‏ خاک Renard, Foster,  (RUSLE)‏ 
(Weesies, McCool, & Yoder, 1997‏ و (Gavrilovic, 1988) *EPM‏ مهم‌ترین این روش‌ها می‌باشند. 

در زمینه کاربرد مدل‌های فرسایش در دنیا تحقیقات متعددی صورت گرفته که در ادامه به برحی اشاره می- 
گردد.(2008) Milevski‏ و )2009( Tazioli‏ به مناسب بودن روش EPM‏ در برآورد فرسایش اشاره نمودند. 
Bahrawi et al.(2016)‏ در بررسی daz‏ فرسایش در عربستان سعودی. تکنیک‌های سنجش‌ازدور و مدل 
RUSLE‏ را تکنیک مناسبی جهت برآورد رسوب اعلام نمودند. Boggs, Devonport, Evans & Puig‏ 
)2001( خطر فرسایش خاک را بر اساس مدل تجدید نظر شده RUSLE‏ و با استفاده از داده‌های مدل 
رقومی ارتفاع" و نقشه واحدهای زمین برآورد نمودند. & Bartsch, Van MiegroetBoettinger,‏ 
Dobrwolski (2002)‏ برای محاسبه فاکتورهای مورد RUSLE 5U‏ جهت تعیین خطر فرسایش خاک کمپ 
cub s‏ از تکنیک‌های GIS‏ استفاده نمودند. بوگیو" و همکاران (۲۰۱۷) فرسایش خاک را در زمان بارندگی- 
های طبیعی در شیب‌های با طول متفاوت بر روی خاک لخت اندازه‌گیری کرده و نتایج آن را با مقادیر 
تخمینی RUSLE ; USLE‏ مقایسه کردند. نتایج نشان داد مقادیر تخمینی بخصوص در مدل USLE‏ کمتر 
از هدر رفت خاک اندازه‌گیری شده بود. 

Moradi, Ghonchepour, «Amiri, Arzani, Farahpour, Chaeichi, & Khajeh Aldin (2009) نیز‎ Ol»! در‎ 
اذعان‎ Ol ul در مناطق خحشک‎ MPSIAC و‎ EPM در مقایسه مدل‎ Nohegar, & Mahmoudi Nejad (2012) 


1- Unversal soil Loss equation 

? - Modified Unversal soil Loss equation 
? - Revised Unversal soil Loss equation 
^ - Erosion Potential Model 

5- DEM 

6 -Bagio et al 
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نمودند. مدل MPSIAC‏ دچار پیش‌برآورد و EPM‏ دچار کم‌برآورد در محاسبه نرخ فرسایش و رسوب ویژه 
می گردد و به ترتیب مدل بدبینانه و خوش‌بینانه قلمداد می‌شوند. )2017( Ansari Lari, & Ansari,‏ و Noori,‏ 
Siadatmousavi, & Mojaradi, (2017)‏ در بررسی مدل EPM‏ مدل فوق را مدل مناسبی برای برآورد 
فرسايش معرفی نمودند. .)2017( Vaezi, Abasi, & Abd Haji Maleki‏ مدل RUSLE‏ همراه ott‏ با 
سنجش‌ازدور و سامانه اطلاعات جغرافیایی را در عرصه‌های زهکش SoS‏ در منطقه نیمهعشک شمال 
غربی ایران ارزیابی کردند. بر اساس نتایج» میانگین هدر رفت خاک برآوردی ۶۳/۷۸ تن در هکتار در سال 
بود که ply VOA‏ بیشتر از میانگین رسوب سالانه AW)‏ تن در هکتار در سال) در عرصه‌های زهکش بود. 
(Mokhtari, Shafiei, & Rahmani, 2019)‏ میزان فرسایش خاک را با استفاده از مدل RUSLE‏ در حوضه 
آبریز نورآباد ممسنی برآورد کردند. نتایج نشان داد میزان خطر فرسایش خاک در سطح حوضه بر حسب تن 
در هکتار در سال از ۸ تا VO‏ متغیر است. ۹270061 Azimi Sardari, Bazrafshan, Panagupulus, & Rafiei‏ 
.)2019( در مطالعه‌ای حال و آینده مقدار فرسایش خاک را در حوضه آبریز سد استقلال میناب با بهره‌گیری 
از مدل RUSLE-3D‏ و سناریوهای تغییر اقلیم ارزیابی کردند. نتایج حاصل از پیش‌بینی بارش در آینده OUS‏ 
دهنده افزايش آن از ۱۶۳/۷ میلی‌متر به ۱۷۸/۱۱ میلی‌متر در دوره آینده و به‌تبع ol‏ افزایش عامل فرسایندگی 
R) abt‏ و درنتیجه افزایش نرخ فرسایش خاک در حوزه سد استقلال میناب است. شرایط سطح خاک مانند 
پوشش گیاهی. زبری و پوشش قطعات سنگی odo‏ مهمی بر روی نرخ نفوذ» فرسایش و رواناب دارد 
(Nouwakpo, Williams, Al- Hamdan, Weltz, Pierson, & Nearing, 2016).‏ 

حوضه آبخیز چم گردلان ایلام نیز یکی از مناطقی است که متحمل فرسایش فزاینده خاک است و بررسی 
فرسایش و بار رسوب در آن به‌عنوان یکی از منابع تأمین آب شرب شهر ایلام از اهمیت بالایی برخوردار 
است. بنابراین» جلوگیری از این فرسایش, به‌عنوان یکی از مهم‌ترین عوامل جهت حفاظت از منابع طبیعی به 
شمار می‌رود )1995 (Morgan,‏ در خصوص مطالعات و تحقیقات انجام شده در حوضه چم‌گردلان ايلام 
می‌توان به رساله کارشناسی ارشد )2001( Azami‏ تحت Ol pe‏ تعیین روش مناسب برآورد بار Glas‏ رسوبی 
در حوضه سد چم‌گردلان ایلام. تخمین فرسایش و رسوب با استفاده از روش‌های AS‏ ژئومورفولوژیکی 
(واحدهای کاری) و EPM‏ و مقایسه آن با آمار خروجی رسوب در حوضه سد چم‌گردلان ایلام Shahbazi‏ 
Asgari, Servati, & Jafari (2008) (1999)‏ تحت عنوان برآورد فرسایش و رسوب حوضه سد gu‏ با 
استفاده از مدل MPSIAC‏ اشاره کرد. 

در دهه‌های اخیر مخاطرات محیطی اعم از طبیعی و انسانی مانند فرسایش خاک ازجمله مسائلی است که در 
سطح جهان و ايران با آن روبرو هستند و سالانه هزینه‌های SUS‏ هم متحمل می‌شوند. ازاین‌ری ابتدا 
مخاطرات محیطی از قبیل فرسایش خاک و تولید رسوب باید شناسایی شوند تا مدیران مربوطه با اطلاعات 
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pls‏ و دقیقی برای مدیریت و JES‏ آن» برنامه‌ریزی مناسبی انجام دهند. ازآنجایی‌که مدل‌های برآورد 
فرسایش و رسوب به‌عنوان ابزاری در راستای شناخت اولویت‌های فرسایش و رسوب حوضه‌های آبخیز 
هستند و مدل‌های EPM‏ و RUSLE‏ بیشترین کاربرد را در مطالعات و تحقیقات آبخیزداری دارد. لذا 
ارزیابی صحت. دقت و نحوه پراکنش طبقات فرسایش و رسوب حاصل از آن‌ها در سطح حوضه می‌تواند 
کمک موثری در نتیجه گیری بهتر از مدل‌ها و به‌کارگیری نتایج آن‌ها در پیشگیری از وقوع مخاطرات زیست- 
محیطی. اولویت‌بندی مناطق برای اقدامات اجرایی و برنامه‌های حفاظتی داشته باشد که در این زمینه مطالعات 
زیادی انجام نشده است. لذا هدف از تحقیق i‏ برآورد فرسایش و میزان رسوب در حوضه آبخیز چم- 
گردلان با استفاده از مدل‌های EPM‏ و RUSLE‏ مقایسه برآورد این مدل‌ها با آمار رسوب مشاهده‌ای: 
بررسی راستی آزمایی و اعتبارسنجی دو مدل نسبت به مقدار واقعی و درنهایت معرفی مدل کارآمد در برآورد 


رسوب در حوضه آبخیز موردمطالعه است. 


منطقه موردمطالعه 

حوضه آبریز چم گردلان در جنوب شرقی شهرستان ایلام در موقعیت بین ۳۰۳ ۱۷۲ PO‏ تا ۳۳ 0۳ £O‏ 
طول شرقی و ۳۲۳ ۲۳ ۳۳٩‏ تا ۵۱۳ YY? £A‏ عرض شمالی قرار دارد. مهم‌ترین آبراهه اصلی در این حوضه 
رود کنجانچم (چم‌گردلان) است که از اتصال دو رود چاویز و ISIS‏ با آبراهه متقاطع تشکیل شده است و 
سد ايلام نیز در محل رودخانه کنجانچم oz)‏ گردلان) احداث شده (Asgari, Servati, & Jafari, Cul‏ 
)2008 سد مخزنی pL!‏ مهم‌ترین منبع تأمین آب شرب شهر ایلام و روستاهای اطراف می‌باشد و بند 
انحرافی امیرآباد در خروجی حوضه و همچنین نیاز به منابع " فراوان در منطقه به علت وجود 
شهرستان مهران و اراضی زراعی مرغوب استان در این منطقه در برنامه‌ریزی توسعه اقتصادی و 
اجتماعی آینده ایس دو شهرستان os‏ بسرایی خواهند داشت. اين حوضه برروی سازندهای 
Asgari, Servati, & Jafari, 2008)‏ 
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شکل اموقعیت کشوری و استانی محد وده موردمطالعه 


۲۱۵۰1۰ National and provincial situation of the study area 


-Y‏ مواد و روش‌ها 

در انجام این تحقیق از مواد زیر با توجه به موضوع. ماهیت تحقیق و روش بکار رفته استفاده شده است: 
الف- آمار و اطلاعات موجود شامل داده‌های هواشناسی» نقشه‌های خاکشناسی» قابلیت اراضی خاک از منطقه 
موردمطالعه در ارتباط با موضوع تحقیق مشتمل بر داده‌ها و اطلاعات مربوط به سازمان هواشناسی و سازمان 
جنگل‌هاء مراتع و آبخیزداری کشور جمع‌آوری گردید. 

ب- DEM‏ حوضه آبخیز با پیکسل سایز ۳۰*۳۰ به‌عنوان نقشه پایه 

ج - نقشه زمین‌شناسی منطقه با مقیاس ۱:۱۰۰۰۰۰ جهت تهیه نقشه سنگ‌شناسی و وضعیت سازند منطقه 

د- تصویر ماهواره‌ای لندست A‏ سنجنده OLI‏ (به تاریخ ۱۳۹۸/۵/۱۵) از سایت Earth explorer‏ جهت 4€ 
تشه cna‏ اتاضتی adl‏ ررد ماله 

»- بسته‌های نرم‌افزاری.5.3 ۳۷ ILWIS3.3‏ و Excel‏ جهت تجزیه‌وتحلیل و استخراج اطلاعات و 
نقشه‌های پایه. 

در این پژوهش از نرم‌افزار ArcGIS10.3‏ جهت ترسیم نقشه‌ها و تحلیل‌های مرتبط به آن؛ ۳۷15.3 جهت 
استخراج و پردازش تصاویر ماهواره‌ای و Excel‏ جهت به دست آوردن روابط رگرسیونی معادلات بهره- 


AXE شده‎ ES 
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روش تحقیق 
RUSLE‏ 
مدل RUSLE‏ یک مدل آماری بر al‏ فرسایش آبی است که با شش فاکتور فرسایشی زیر رابطه دارد 


:(Wischmeier & Smith, 1978) 
A=R.K.L.S.C.P A) 

در این رابطه» A‏ میانگین فرسایش خاک در واحد سطح (تن در هکتار در R (Jle‏ فاکتور فرسایندگی باران 
(مگاژول میلی‌متر در هکتار ساعت K (UU.‏ فاکتور فرسایش‌پذیری خاک (تن ساعت بر مگاژول میلی‌متر) 
فاکتور طول شیب. S‏ فاکتور درجه شیب. C‏ فاکتور عملیات مدیریتی و پوشش و P hle‏ فاکتور عملیات 
حفاظتی می‌باشد. مقادیر CIS L‏ و P.‏ بدون واحد می‌باشند. در این پژوهش. روش بررسی fale IS‏ استفاده 
از مدل RUSLE‏ در GIS bes‏ می‌باشد. پاراگراف‌های بعدی نیز درباره نحوه برآورد فاکتورهای PCKR‏ 
و LS‏ که به ترتیب از داده‌های بارش نقشه‌های خاک پردازش تصاویر ماهواره‌ای (فاکتورهای € و P‏ و 


فاکتور فرسایندگی باران (R)‏ 

فرسایندگی باران بە‌صورت قدرت تراکمی باران در بروز فرسایش تعریف می‌شود )1990 (Lal,‏ معمول‌ترین 
شاخص فرسایندگی hb‏ فاکتور R‏ مربوط به مدل‌های (Yu & Rosewell, 1996) USLE‏ و RUSLE‏ 
(Renard & Freimund, 1994)‏ می‌باشد. بررسی مقالات علمی نشان می‌دهد. در بسیاری از نقاط di»‏ 
فاکتور *آهمبستگی بالایی با فرسایش خاک دارد. فاکتور R‏ برای دوره‌های مختلف از حاصل‌ضرب انرژی 
جنبشی باران (B)‏ در حداکثر شدت بارش ۳۰ دقیقه‌ای (130) به دست می‌آید. ازآنجایی که نمودار بارندگی و 
داده‌های تفصیلی رگبار (شدت بارندگی) به‌ندرت در ایستگاه‌های هواشناسی موجود می‌باشند. غالباً از مقادیر 
متوسط بارندگی سالیانه و ماهیانه» جهت برآورد فاکتور R‏ در مدل‌های USLE‏ و RUSLE‏ استفاده می‌شود. 
برای محاسبه فاکتور R‏ بعد از تعیین ایستگاه‌های شاخحص در منطقه موردمطالعه» بارندگی ماهیانه و سالیانه 
در این ایستگاه‌ها و در دوره زمانی موردمطالعه بازسازی گردید (بر اساس اطلاعات ایستگاه ايلام که 
کامل‌ترین و بیشترین آمار را داشته است). در مرحله بعد با استفاده از معادله ذیل» شاحص فورنیر و فاکتور R‏ 


برای تمام ایستگاه‌ها به دست آمد. رابطه شاحص Fog‏ به‌صورت eal d)‏ 


12 p;2 
154 Pi 
F = کن‎ 


در این معادلی pi‏ متو سط بارندگی (میلی‌متر) در ماه 1و 0 متوسط بارند گی all.‏ (میلی‌متر) AC aol‏ در این 
مطالعه. با استفاده از معادله Y)‏ شاخحص فورینر برای تمامی ایستگاه‌هاه محاسبه گردید و سپس با جایگذاری 


۳۴۶ نشریه جغرافیا و مخاطرات محیطیء جلد Y‏ شماره ۲ تابستان ۱۶۰۳ 


شاخحص فورنیر (رابطه (Y‏ در روابط زیر که توسط )1994( Renard & Freimund‏ برای مناطق فاقد داده‌های 
تفصیلی رگبار (شدت بارندگی) پيشنهاد شده. مقدار فاکتور R‏ برای ایستگاه‌های شاخحص برآورد شد. 
اگر R — factor = (0.07397 x 2 (۳) F«55mm‏ 
اگر R — factor = (95.77 — 6.081 x F + 0.4770 x F?/17.2 (£) F>=55mm‏ 


فاکتور فرسایش پذیری خاک BK)‏ 

فاکتور فرسایش‌پذیری خاک (K)‏ سرعت فرسایش‌پذیری خاک به ازای واحد شاخص فرسایش OL‏ است 
که در GL‏ استاندارد اندازه‌گیری می‌شود و اغلب با استفاده از خحصوصیات خاک تعیین می‌شود. ax S|‏ 
فاکتور UK‏ بافت خاک» میزان ماده آلی SE‏ نفوذپذیری و سایر فاکتورها رابطه دارده اما عمدتاً با نوع تیپ 
SE‏ همبستگی بالایی دارد. فاکتور فرسایش‌پذیری SE‏ در مدل RUSLE‏ از طریق نوموگراف تعیین 
می‌شود. این نوموگراف متشکل از پنج پارامتر پروفیل خاک می‌باشد. رابطه جبری این نوموگراف به‌صورت 

K = 2.73 x 10-5m*14(12 — a) + 3.25 x 1072 (b — 2) + 2.5 x 107-266 — 3) 

در این معادله» M‏ قطر o>‏ (درصدسیلت + درصد شن خیلی ریز)»(درصد رس - acd\es‏ درصد ماده b cdl‏ 


کد ساختمان خاک و 6 کلاسه نفوذپذیری پروفیل خاک می‌باشند. در این مطالعه تمام پارامترهای اشاره شده 


فاکتورهای طول شیب (D)‏ و درجه شیب (S)‏ 

S 4L clay st‏ در مدل RUSLE‏ تأثیر توپوگرافی را بر روی فرسایش خاک نشان می‌دهد. Col‏ شده که 
افزايش در طول شیب و درجه شیب می‌تواند سرعت جریان آب را در سطح زمین افزايش دهد و باعث 
تشدید فرسایش خاک شود )1994 (Haan, Barfield, & Hayes,‏ و از اینها گذشته فرسایش خالص خاک 
نسبت به تغییرات در درجه شیب در مقابل طول شیب حساس‌تر (McCool, Brown, & Foster, cul‏ 
(1987. طول شیب فاصله بین نقطه بالایی آغاز سراشیبی تا نقطه‌ای است که شیب به حداقل خود می‌رسد و 
در این نقطه رسوبات ته‌نشین می‌شوند. تأثیرات ویژه توپوگرافی روی فرسایش خاک توسط فاکتور LS‏ 
(بدون واحد)(طول شیب (L)‏ و درجه شیب (S)‏ برآورد می‌شود. برای محاسبه فاکتورهای بآ و S‏ از 
فرمول‌های زیر استفاده می‌شود(1993 (Moore & Wilson,‏ 

s < [ey (v) L = 1.4-5 CU 

در این معادله As‏ مساحت aul‏ بالادست در واحد عرض خطوط توپوگرافی (m/m)‏ و B‏ زاویه شیب 
برحسب درجه است. در این مطالعه. نقشه فاکتور 18 از DEM «zz‏ با قدرت تفکیک ۲۰ متر از وب سایت 
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سازمان نقشه‌برداری امریکا به آدرس http://earthexplorer.usgs.gov‏ دانلود و با استفاده از اکستنشن 
0۴۲ آرک‌جی‌آی اس ۱۰.۳ استخراج شد که براساس کار Moore & Wilson(1993)‏ توسعه ail,‏ 
است. 
فاکتور مدیریت پوشش گیاهی (C)‏ 
فاکتور مدیریت پوشش گیاهی (C)‏ نشان‌دهنده تأثیر فعالیت‌های کاشتی در مدیریت کشاورزی و اثر 
پوشش‌های درختی» درختچه‌ای, علفی و زمینی روی کاهش فرسایش خاک می‌باشد. با افزایش پوشش 
گیاهی. فرسایش خاک کاهش می‌یابد. )1994( Biesemans, Meirvenne, &  Benkobi, Trlica, & Smith‏ 
Gabriels (2000)‏ گزارش کردند. فاکتور پوشش گیاهی به همراه درجه شیب و طول شیب در فرسایش 
خاک و تولید رسوب حساسیت بیشتری دارند. در مدل RUSLE‏ فاکتور پوشش گیاهی (C)‏ معمولاً بر 
اساس معادلات تجربی تعیین می‌گردد. اما پرکاربردترین معیار رویش پوشش گیاهی. شاخص تفاضل پوشش 
گیاهی نرمال شده است که با استفاده از فناوری سنجش ازدور به دست می‌آید. برای تضویر لندست این 
شاخحص به‌صورت bd‏ به دست می‌آید: 
)^( 


NIR-IR 
NIR+IR 


این شاخص معرف انعکاس انرژی خورشیدی از سطح زمین است که انواع شرایط پوشش گیاهی را نشان 
می‌دهد. NDVI polie‏ بین ۱- و ۱+ در نوسان است. زمانی که پاسخ طیفی اندازه گیری‌شده از سطح زمین 
Gly‏ هی کو ales Je ath‏ باک NDVI pot‏ په ضقن رھک می قود پر gal.‏ مال Bile)‏ 
فعالیت فتوسنتزی) در بخش طیفی مادون‌قرمز نزدیک (NIR)‏ (باند ۵ لندست) در مقایسه با بخش طیف 


NDVI = 


خواهد شد. مناطق با پوشش گیاهی کم یا بدون پوشش گیاهی از قبیل مناطق شهری و اراضی yh‏ معمولا 

NDVI polis‏ بین ۰/۱+ و ۰/۱- را نشان می‌دهند. Lay!‏ و منابع آبی مقادیر منفی یا صفر را OUS‏ می‌دهند. 

شاخص تفاضل پوشش گیاهی نرمال شده (NDVI)‏ برای محاسبه داده‌های طیفی زمینی استفاده شده و نتایج 

él‏ نشان داده که این شاحعص. همبستگی YL‏ با زیست‌توده سطح زمین دارد. بعد از انجام تبدیل حطی 
C = (1— NDVD/2 (4)‏ 


درنتیجه مقدار € در هر سلول می تواند موردمحاسبه قرار گیرد. 


۱۶۰۳ تابستان‎ Y جغرافیا و مخاطرات محیطی» جلد ۰۱۳ شماره‎ a iS YFA 


فاکتور عملیات حفاظتی (P)‏ 
در این مطالعه» poles‏ فاکتور P‏ از طریق طبقه‌بندی مجدد هر تیپ پوشش زمین و بر اساس 
oslo polis‏ شده در جدول‌های Y‏ و Y‏ تعیین شدند )2006 (Deore,‏ با توجه به اینکه در حوضه 
تحت بررسی هیچ‌گونه عملیات اصلاحی و مدیریتی صورت نگرفته. ارزش P‏ تنها براساس کاربری و پوشش 
اراضی مختلف تعیین شده است (جدول ۱ و (Y‏ 

حدول ۱- ارزش P‏ برای کاربری و پوشش اراضی مختلف )2006 (Deore,‏ 


Table 1- P value for different land use and coverage 


کلاس کاربری/پوشش اراضی ارزش عددی P‏ 
اراضی ak‏ ۱ 
نیشکر IW‏ 
گندم n‏ 
جنگل تنک IA‏ 
اراضی آیش ۱ 
جنگل با تراکم متوسط IA‏ 
جنگل باز IA‏ 
بستر رودخانه ۱ 


(Deore, 2006) برای عملیات حفاظتی مختلف‎ P ارزش‎ -۲ Jj 


Table 2- P value for different protective measurement 


شیب ارزش عددی P‏ 

(درصد) کشت روی byl‏ تراز کشت نواری PT‏ بدون اقدامات حفاظتی 
E ۰/۸ "T‏ = \ 

\ -= UN ۰/۹ ۱-۲ 

۱ A ۰/۵ ۰/۵ y-V 

\ ARI VAN ۰/1 ۷-1۲ 

\ /M ZH ۰/۸ \Y-\A 

\ /M ۰:۵ ۰/۹ 1۸-٤ 


محاسبه نسبت تحویل رسوب (SDR)‏ 
در اين مطالعه سه روش بر اساس مساحت (0854,1975 1975 (Vanoni, 1975 Boyce,‏ به‌منظور 
محاسبه نقشه نسبت تحویل رسوب مورد استفاده قرار گرفتند. معادلات مربوط به این سه روش به شرح زیر 
۱- مدل وانونی (ارائه شده در سال (VAVO‏ به شرح رابطه ۱۰ 

SDR = 0.47244 925  (y.) 
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که A‏ = مساحت حوضه (کیلومترمربع) 
۲- مدل ارائه شده USDA bu y‏ در سال ۲ به شرح رابط ۱۱ 
SDR = 0.5656۸ 911. (1)‏ که ۸ = مساحت حوضه (کیلومترمربع) 
۳- مدل بویس (VAVO)‏ که به شرح رابطه ۱۲ است 
SDR = 0.37508: 92382 . (v)‏ که A‏ = مساحت حوضه (کیلومترمربع) 


مدل EPM‏ 
در روش EPM‏ چهار مشخصه شامل ضریب فرسایش حوضه آبخیز QU)‏ ضریب استفاده از زین (Xa)‏ 
ضریب حساسیت سنگ و خاک به فرسایش (Y)‏ و شیب متوسط حوضه (D‏ در واحدهای همگن 
موردبررسی قرار گرفت. در هر یک از واحدهای همگن بر اساس این چهار «le‏ مقدار (Z)‏ یا ضریب 
شدت فرسایش از رابطه ۱۳ محاسبه گردید )1988 (Gavrilovic,‏ 

Z = ۷.۵) ۱ ۲۳ که در آن‎ (Y) 
عامل حساسیت خاک و سنگ به فرسایش‎ : Y 
عامل استفاده از زمین‎ : Xa 
عامل فرسایش‎ : ۴ 
عامل شیب زمین به درصد‎ :1 
(Gavrilovic, 1988) (Y در جدول استاندارد (جدول‎ ol با محاسبه ضریب شدت فرسایش و قرار دادن مقدار‎ 
تنظیم شده است» طبقه فرسایشی هر یک از‎ (Z) که بر اساس شدت فرسایش و ضریب شدت فرسایش‎ 
واحدهای همگن مشخص می‌شود.‎ 

جدول ۳- طبقه‌بندی شدت فرسایش 


Table 3- Erosion intensity classification 


کلاس فرسایش شدت فرسایش محدوده Z‏ مقدار متوسط Z‏ 
M‏ خیلی شدید )> VYo‏ 
IV‏ شدید ۵ >\ ۰/۸0۵ 
11 متوسط ۱ < ۰/۷ 00/* 
I‏ کم ۲ < ۰/1 Zu‏ 


1 خیلی کم Z<+/\4‏ ۰/۱ 


۱۶۰۳ شماره ۲ تابستان‎ Y نشریه جغرافیا و مخاطرات محیطیء جلد‎ Yo. 


در روش EPM‏ برای تخمین متوسط سالانه فرسایش ویژه که در حقیقت میزان فرسایش را در واحد سطح 
(هکتار یا کیلومترمربع) حوضه را نشان می‌دهد. از رابطه ۱۶ استفاده می‌شود. 
Wsp = T.H.m.Z15 13‏ 

که در آن: 
H‏ بارندگی متوسط سالانه به میلی‌متر» T‏ عدد ثابت پی برابر با ۳/۱۶ 
T‏ ضریب درجه حرارت که از رابطه ۱۵ به‌دست می‌آید )1988 (Gavrilovic,‏ 

T = (Š + 0.1)°5 (\0) 

که در ol‏ 

= میانگین درجه حرارت سالانه در حوضه به درجه سانتی‌گراد می‌باشد. بعد از محاسبه فرسایش ویژه مقدار 
فرسایش سالانه برای حوضه موردمطالعه با استفاده از رابطه VV‏ محاسبه می‌شود. 

Ws = Wsp X A ("9‏ 
که در آن 
و : فرسایش کل حوضه بر حسب مترمکعب در سال» 
۸ مساحت حوضه بر حسب کیلومترمربع 
Wsp‏ محاسبه شده در رابطه فوق بیانگر مقدار SLE‏ است که از بستر خود جدا کرده و انتقال یافته است؛ اما 
همه این مواد فرسایش EL‏ به خروجی حوضه نمی‌رسد و آن مقدار خاک فرسایش al‏ که به رسوب 
خروجی تبدیل خواهد شد. بستگی به متغیرهای مختلفی دارد که در واقع بر نسبت رسوب‌دهی حوضه jd‏ 
می‌گذارند. در این مدل از ضریبی برای تبدیل میزان فرسایش به رسوب استفاده می‌شود که ضریب رسوب- 


T‏ حوضه نامیده می‌شود. 


محاسبه ضریب رسوبدهی حوضه 
با اعمال این ضریب معین می‌شود که چه ميزان مواد رسوبی از نقطه خروجی حوضه خارج خواهد شد. 
ضریب رسوب‌دهی در روش EPM‏ براساس (Gavrilovic, 1988) ۱۷ 4b;‏ و با استفاده از عواملی مثل 
محیط و ارتفاع و طول حوضه به دست می‌آید: 

Ru = 4 x (Ox D)?5/[L + 10] Qv) 
که در این رابطه:‎ 
اختلاف‎ D کیلومتر‎ Goel طول حوضه‎ L ضریب رسوب‌دهی حوضه» 0: محیط حوضه به کیلومتر:‎ Ru 
ارتفاع بین ارتفاع متوسط و ارتفاع نقطه خروجی حوضه به کیلومتر که از رابطه ۱۸ محاسبه می‌شود:‎ 
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D = Day — Do(\A) 
که در این رابطه:‎ 
ارتفاع متوسط حوضه آبخیز بر حسب کیلومتر‎ Day 
را پس از محاسبه» در رابطه‎ D ارتفاع نقطه خروجی در رودخانه بر حسب کیلومتر است؛ بنابراین مقدار‎ Dy 
محاسبه می‌شود‎ ۱٩ قرار داده و ضریب رسوب‌دهی حوضه محاسبه می‌شود. سپس رسوب ویژه از رابطه‎ ۵ 
.(Gavrilovic, 1988) 
Gsp 2 Wsp x Ru (44) 
که در این‎ 
تولید رسوب ویژه» بر حسب مترمکعب در کیلومترمربع در سال‎ : 9 
فرسایش ویژه. بر حسب مترمکعب در کیلومترمربع در سال‎ : 0 
(Gavrilovic, 1988) ضریب رسوب‌دهی حوضه که از رابطه ۱۷ به دست می‌آید‎ : Ru 
و درنهایت رسوب کل حوضه از رابطه ۲۰ محاسبه می‌شود:‎ 
و و6۵‎ (Ys) 
که در این رابطه:‎ 
رسوب کل بر حسب مترمکعب در سال در حوضه‎ : Gs 
مساحت حوضه آبخیز بر حسب کیلومترمربع می‌باشد.‎ A 
و‎ EPM در این تحقیق برای برآورد شدت فرسایش و پتانسیل رسوبدهی حوضه آبخیز موردمطالعه از روش‎ 
و زمین مرجع کردن‎ ArcGIS استفاده شده است. بعد از وارد کردن نقشه‌ها به محیط نرم‌افزار‎ GIS تکنیک‎ 
EPM رقومی و پایگاه اطلاعاتی برای آن‌ها ایجاد شد. سپس با استفاده از مدل‎ GLE لایه‌های مورد‎ daol 
عوامل مؤثر در فرسایش در واحدهای همگن (شکل ۲) ارزیابی و امتیازبندی شد. در این روش که در‎ 
L ضریب استفاده از زمین‎ QU) حقیقت ضریب شدت فرسایش است. چهار عامل ضریب فرسایش حوضه‎ 
در‎ (D و شیب متوسط حوضه آبخیز‎ (Y) ضریب حساسیت سنگ و خاک به فرسایش‎ (Ka) کاربری اراضی‎ 
واحدهای همگن و درنهایت در کل حوضه موردمطالعه قرار گرفت. در هریک از واحدهای همگن, بر اساس‎ 
رسوب‎ JS یا ضریب شدت فرسایش. متوسط سالانه فرسایش ویژه» فرسایش‎ Z چهار عامل فوق. مقدار‎ 
ویژه و رسوب کل تولید شده در طول سال در سطح واحدهای همگن و بالاخره در کل حوضه آبخیز محاسبه‎ 
گردیده است.‎ 
همان‌طور که ذکر شد. برای سهولت و دقت کار در به دست آوردن میزان فرسایش و تولید رسوب. ابتدا کل‎ 


حوضه آبخیز چم گردلان به واحدهای همگن تقسیم شد (۱۵ واحد همگن). سپس عوامل ja‏ در تولید 


۳۵۲ نشریه جغرافیا و مخاطرات «dama‏ جلد Y‏ شماره ۲ تابستان ۱۶۰۳ 


رسوب در هر یک از واحدهای همگن موردمطالعه قرار گرفت. در مرحله Gel‏ میزان فرسایش و تولید 
رسوب با استفاده از فرمول‌ها و ضرایب مربوطه به دست آمد که در زیر به بررسی هر یک از عوامل چهارگانه 


خواهیم پرداخحت. 
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شکل -Y‏ واحدهای همگن 


Fig.2. Homogeneous units 


ارزیابی کارایی مدل‌ها 
cis Lol‏ مدل ely‏ اطمیتان از lasgul yy Gone‏ ضروری است» ارژبانی مذل‌ها با اسفاده از داده‌های 
برآوردی رسوب و داده‌های حاصل از اندازه‌گیری رسوب در خروجی حوضه انجام گرفت. همچنین کارایی 
مدل که اغلب برای کالیبراسیون صحت سنجی و برآورد دقت مدل‌ها مورد استفاده قرار می‌گیرد از رابطه زیر 
تعیین شد (رابطه ۲۱). پس از برآورد رسوبات توسط مدل مذکور برای ارزیابی کارایی مدل بکار رفته از 
رابطه زیر استفاده گردید. 

Rg = (ISo — Sg|/S5) x 100 (3)‏ 
2=RE‏ صد خطای نسبی «JJ‏ 
=So‏ مقدار رسوب مشاهده‌ای, 
=SE‏ مقدار رسوب برآوردی 
cual essa Ls 2 Re dd‏ کش City‏ بای aes Sle atl‏ کل که مان 
متغیر موردنظر را کم یا زیاد برآورد می کند. از نظر تئوری هرگاه این دو مقدار برابر با صفر شوند. نمایانگر 
این مطلب است که دقت روش صد در صد بوده و مقدار برآورد شده یک کمیت» دقیقاً برابر با مقدار واقعی 


آن است. 
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نتایج و بحث 

RUSLE Ja 

شاخص‌های فورنیر و مقادیر فرسایندگی باران با استفاده از معادلات Q0‏ (۳) و E)‏ برای ۱۸ ایستگاه در طی 
دوره آماری ۲۲ ساله (۱۳۹۸-۱۳۷) محاسبه شد. polis‏ محاسبه شده شاحص phy gd‏ و فاکتور R‏ در جدول 
OUS 6‏ داده شده است. ابتدا بر اساس رابطه بین ارتفاع و بارش و سپس بر اساس رابطه بارش و فاکتور R‏ 
نقشه فاکتور ۸ حوضه چم‌گردلان در محیط نرم‌افزار ILWIS3.3‏ به دست آمده است (شکل (Y‏ مقدار 
متوسط فاکتور GIR‏ ۱۳۱ تا ۳۹۲ (مگاژول میلی‌متر در هکتار ساعت سال) متغیر بوده و مقدار متوسط و 
انحراف معیار آن به ترتیب ۲۱۶ و ££ (مگاژول میلی‌متر در هکتار ساعت سال) می‌باشد. فرسایندگی باران در 
شمال‌شرقی و جنوب‌غربی حوضه موردمطالعه که ارتفاع بالاتری داشتند بیش‌تر بوده که با رنگ قرمز نشان 
داده شده است L‏ به‌عبارت‌دیگر. کاهش ارتفاع و باران رابطه مستقیمی با کاهش فاکتور R‏ دارد و از شمال- 


غرب و جنوب‌غرب به مرکز حوضه داشته است. 


فاکتور فرسایشپذیری خاک K)‏ 
جهت محاسبه فاکتور K‏ اطلاعات مورد نیاز از نقشه‌های تفصیلی حوضه چم گردلان و لایه‌های زمین‌شناسی؛ 
خاک‌شناسی استفاده شده که متوسط مقدار K‏ در حوزه موردمطالعه بین ۰/۲۸ تا ۰/۹۱ متغیر بوده میانگین و 
انحراف معیار Ol‏ به ترتیب ۰/۳۵ و ۰/۰۵ تن ساعت بر مگاژول میلی‌متر می‌باشد. شکل ٤‏ توزیع مکانی 
elas‏ تا برغ اک رز ان das‏ 

جدول -٤‏ محاسبه و برآورد ۴ و 8 برای ایستگاه‌های بارندگی 


Table 4- Calculation and estimation of F and R for rainfall stations 


R F نام ایستگاه طول عرض ارتفاع  بارندگی‎ 
بے (میلی‌مت) مگاژول میلی‌متر در هکتار ساعت‎ delle پان‎ 
سال‎ (UTM) (UTM) 
YEY ۱۱3/۵۶ 1۳/۵۰ ۰ 40۸14 MAEVA 3583! 
\AV/00 AVAY ۳۱۳۸۷ WS Yve £YYA When صالح آباد‎ 
YA0/£* AVIAN £A\ qv» ۱۳۷۱۷۰ ۳۷۹ سراب کلان‎ 
۱/۸ ۸۱/۳۲۱ ۷۷/۰ ۰ ۱۳۹۹۲۸۰ 1014۰۰ ds 
طالقانی‎ 
۱۲۱۳۱ AWAIT ۰ ۵۲۲/۵۰ ۰ TVITAAT 1764 طولاب‎ 
\AVAA ۸۷/۸ — 1۹۹/۵۰ 140 ۱۳۳۹۹۰۹ 1011 قلندرعلیا‎ 


9 1۱۲ 0۹/7۰ ۱۳۵۰ TVA KAUSE? میشخاص‎ 


۱۶۰۳ شماره ۲. تابستان‎ VY جغرافیا و مخاطرات محیطی. جلد‎ 4 pia Yor 


was ۸0/۸٦ 4٤ ۱۳۸۰ Wry MAT: کارزان‎ 
YAY'/00 AWE tvé NV Y" £04 MYYSY گناباد‎ 
۳/۸ T/A ۰ ۰ rmt EVEL درگه‎ 
۹710 ۸۹/۵ — EVVIV ۸۰ YY) 0008 WAY قجر‎ 
VAAT ۵۸/۱۸ ۲۸۸/۷۰ ۰ ۰ Y11444) woo بان رحمان‎ 
VE ۹۳۸۷/۲ 1۸۹ ME ۱۳۰ ۳۹۹۰ سیاه آب‎ 
VIE ۱۱/۹۰ 0۹۷ ۱۱۶۰ RAVA ۳۷/۳ گل گل‎ 
۱۳/۳۰ VVA\ TEKI: ow ۱۳۶ IVE! گلان‎ 
۳1/00 NIVEL AVIV ۷ ۰۸۹ T1 ايلام‎ 
AVEN MY ۲۳۱ Ye 11017 "AV مهران‎ 
۱۳۸۳ VT 011/0 14 PVs sas TYE اما‎ 
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شکل -Y‏ نقشه عامل فاکتور R‏ 


Fig.3. R factor factor map 


K) SE نقشه فاکتور فرسایش پذیری‎ -٤ شکل‎ 
Fig.4.Map of Soil erodibility factor (K) 


فاکتورهای طول شیب (D)‏ و درجه شیب (S)‏ 
نقشه توپوگرافی منطقه موردمطالعه به‌منظور ایجاد طول شیب و درجه شیب (LS)‏ مورد استفاده قرار گرفته 
AR‏ فاکتور LS‏ به‌وسیله معادله‌های مربوطه با استفاده از نقشه DEM‏ حوضه و با در نظر گرفتن کنش 


متقابل بین توپوگرافی و تجمع جریان آب محاسبه شد. همان‌طور که از شکل ٦‏ برمی‌آید. فاکتور LS‏ در 
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منطقه موردمطالعه بین ۰ تا ۱۶۷ در سطح پیکسل متغیر بوده و مقدار متوسط آن ۲ می‌باشد. مقادیر LS‏ در 


PON ا ا روات ها و الها ی تر‎ E ا ا‎ NEPOS 


فاکتور مدیریت پوشش (C)‏ 


پوشش متراکم وجود دارد که با رنگ روشن تر نشان داده شده است. 
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(LS) شکل ۵ نقشه فاکتور طول شیب و درجه شیب‎ 
۲۱2۰5۰ Map of Slope length factor and slope degree (LS) 


Fig.6. Coverage Management Map, Invoice (C) 


- فاکتور عملیات حفاظتی (P)‏ 
با توجه به اینکه در حوضه تحت بررسی هیچ‌گونه عملیات اصلاحی و مدیریتی صورت نگرفته؛ فاکتور 


حدول ۵ ارائه شده است. 
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(P) شکل ۷- نقشه عملیات حفاظتی فاکتور‎ 
۲۱۵۰7۰ Map of factor protection operations (P) 
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جدول ۵- آمار مربوط به فاکتورهای معادله RUSLE‏ در حوضه چم گردلان ایلام 


Table 5- Statistics related to the factors of the RUSLE equation in the Cham- 
Gardlan basin of Ilam 


P فاکتور‎ C فاکتور‎ LS فاکتور‎ K فاکتور‎ R فاکتور‎ 
\ ۳۵ ۰ ۰/۳۸ ۱۳۹ حداقل‎ 
۱ De ۱1۷ Za ۳۹3 ES 
\ ۰/۵۱ Y ۳۵ u میانگین‎ 
' ۳ ۹ TID £t انحراف معیار‎ 


فرسایش سالانه خاک 

متوسط فرسایش GUL‏ خاک به‌وسیله ضرب کردن فاکتورهای فرسایندگی (R)‏ فاکتور فرسایش‌پذیری AK)‏ 
فاکتور توپوگرافی (LS)‏ فاکتور مدیریت پوشش گیاهی (C)‏ و فاکتور عملیات حفاظتی (P)‏ در محیط 
نرم‌افزاری 10.3 ArcGIS‏ محاسبه شد. نقشه به دست آمده از این مدل» در شکل ۸ نشان داده شده است. 
polis‏ فرسایش خاک در منطقه موردمطالعه بین صفر و ۶۷۳۰ تن در هکتار در سال تعیین شد. در مرحله بعد 
نقشه فرسایش سالانه خاک به ۵ کلاس خطر طبقه‌بندی شده است polis 4 JSS)‏ فرسایش خاک در منطقه 
موردمطالعه بین صفر Og)‏ فرسایش) تا 1۷۳۰ تن در هکتار در سال در سطح پیکسل متغیر می‌باشد. متوسط 
و انحراف معیار OF‏ نیز به ترتیب برابر با ۸7/۵۰ و ۳۶۳/۰۰ تن در هکتار در سال تعیین شد. نقشه به دست 
آمده از این مدل. در شکل ۸ نشان داده شده است. در مرحله بعد نقشه فرسایش سالانه خاک به ۵ کلاس 
خطر طبقه‌بندی شده است (شکل 4( )2008 (Yue-Qing, Shao, & Kong,‏ فراوانی هر یک از کلاس‌های 
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خطر در شکل oU‏ داده شده Aral‏ بیشترین قسمت حوضه موردمطالعه. جزء طبقه فرسایشی خیلی کم 
تا کم V0)‏ درصد) قرار داشت که عدا این کلاس‌ها در بخش میانی حوضه واقع شده بودند. حدود ۲۰ 
درصد حوضه نیز تحت خطر زیاد تا خیلی زیاد قرار داشت که ازاین‌بین» بیشترین مناطق تحت خطر فرسایش 
پژوهش همچنین OUS‏ داد. فاکتور LS‏ با ضریب همبستگی ۱ بیشترین تأثیر را در برآورد فرسایش سالانه 
خاک توسط مدل RUSLE‏ داشته است (شکل OY‏ این در حالی است که همبستگی بسیار پایینی بین سایر 
فاکتورها و فرسایش SUL.‏ خاک وجود داشت (جدول CV‏ 
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شکل ۸- نقشه متوسط سالانه فرسایش 


Fig.8. Annual average erosion map 


شکل -٩‏ نقشه کلاس فرسایش خاک 


۲12.9. Soil erosion class map 
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شکل ۱۱- رابطه فاکتور LS‏ با متوسط سالانه فرسایش خاک شکل ۱۰- فراوانی کلاس‌های خطر فرسایش 
Fig.10. Frequency of erosion risk classes‏ 


Fig.11. Relationship between LS factor and average annual soil erosion 


۳2۸ نشریه جغرافیا و مخاطرات محیطیء جلد Y‏ شماره Y‏ تابستان ۱۶۰۳ 


حدول M‏ — میزان همبستگی و روابط رگرسیونی بین مقدار فرسایش سالانه خاک و پارامترهای مدل RUSLE‏ 


Table 6- 
RUSLE سالانه خاک فاکتورهای‎ cub 
R r = 0.0024, y = -0.3526x + 195.16 
K r= 0.0151, y = 270.91x + 198.17 
LS 1= 0.8111, y = 20.337x + 45.145 
C r= 0.0866, y = 503.16x - 67.268 
P T= 0.0993, y = 1642.6x - 1071.3 


بار رسوب 

در این مطالعه. متوسط polis‏ نسبت تحویل رسوب بین ۰/۱۲ تا ۰/۳۶ و متوسط بار رسوب در خروجی 
حوضه بین ۵/۳۱ تا ۱۷/۳۷ تن در هکتار در سال متغیر بود (جدول (V‏ 

میزان رسوب مشاهده و محاسبه شده بر اساس VV‏ سال آمار ایستگاه هیدرومتری با استفاده از منحنی‌های 
سنجه (حد وسط دسته‌ها) به‌طور متوسط Sales‏ ۱۱/۵۸ تن در هکتار در سال )2001 (Azami,‏ بوده و مقدار 
رسوب برآورد شده (متوسط رسوب برآورد شده) بر اساس مدل RUSLE‏ با استفاده از روش نسبت تحویل 
رسوب Joles USDA‏ ۱7/۳۷ تن در هکتار در سال می‌باشد. در ضمن مقدار بار رسوبی برآورد شده در این 
حوضه توسط )1999( Shahbazi‏ با استفاده از روش EPM‏ برابر ۱۸/۹۰ تن در هکتار در سال (آمار مشاهده- 
ای معادل VV‏ تن در هکتار در سال) و(2008) L Asgari, Servati, & Jafari,‏ استفاده از مدل MPSIAC‏ 
برابر ۱۶/۹۸ تن در هکتار در سال می‌باشد (آمار مشاهده‌ای ۱۱/۵۸ تن در هکتار در سال). بنابراین اختلاف 
بین مقدار رسوب مشاهده شده با مقدار برآورد شده توسط مدل RUSLE‏ در این مطالعه برابر ۰/۲۱ تن در 
هکتار در سال است. به‌عبارت‌دیگر مقدار رسوب برآورد شده حدود ۱ در صد کمتر از مقدار رسوب اندازه 
گیری شده (مشاهده شده) (براساس محاسبه خحطای نسبی) می‌باشد. 

بنابراین اختلاف معنی‌داری بین نتایج این کار و آمار واقعی مشاهده نمی‌شود؛ بنابراین ازآنجایی که مقدار بار 
رسوب برآوردی با استفاده از روش USDA‏ (جدول V‏ تقریباً نزدیک به مقدار بار رسوب ثبت‌شده در 
ایستگاه رسوب‌سنجی حوضه چم‌گردلان ایلام VVOA)‏ تن در هکتار در سال) بود )2001 (Azami,‏ ازاین‌ری 
این روش به‌عنوان بهترین روش برای محاسبه نسبت تحویل رسوب انتخاب گردید و براساس Ol‏ نقشه 
نسبت تحویل رسوب و بار رسوب تهیه شد (شکل‌های ۱۲ و ۱۳). در این روش» متوسط نسبت تحویل 


رسوب و رسوب سالانه از ۵۹۵ ۲ و ۰ تا 980۸ تن در سال به ترتیب متغیر است. 


سال سیزدهم آرخی و همکاران. مقایسه مدل‌های RUSLE‏ و EPM‏ در برآورد فرسایش و تولید رسوب با استفاده از ... YO‏ 


جدول ۷- روش‌های نسبت تحویل بار رسوب. مقادیر نسبت تحویل و بار رسوب در حوزه موردمطالعه 


Table 7- Sediment load delivery ratio methods, delivery ratio values and sediment load in 
the study area 


روش‌های تحویل بار رسوب Ol jee‏ نسبت تحویل بار رسوب (تن در هکتار در سال) 
وانونی Um‏ ۱۱/۳۲ 
Zu USDA‏ ۱/۳۷ 
Ww w Mee Mee 11 w‏ 
N w Were Levee v 1" Ww‏ 
Eo d‏ £ 
È‏ | | ن مرز منطقه مورد مطالعه 3 
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Fig.12. Sediment delivery ratio map 


Fig.13. Sediment load map 


محاسبه فرسایش و رسوب با مدل EPM‏ 

در حوضه موردنظر, برای تعیین میزان فرسایش خاک و همچنین تولید رسوب از مدل EPM‏ استفاده گردید. 
با توجه به نتایج به دست آمده از مطالعات پایه و بازدیدهای میدانی انجام شده. ضرایب تعدادی از 
پارامترهای مدل برای منطقه موردمطالعه به شرح جدول ۸ می‌باشد. در این قسمت. میانگین چهار ضریب؛ 
موردبررسی و تجزیه‌وتحلیل قرار می گیرد. نتایج چهار ضریب در فرسایش خاک در جدول ۸ مشخص شده 
است و سپس این چهار M‏ (شکل OE‏ ۸۱۵ ۱۱ و OV‏ در محیط نرم‌افزار ArcGIS10.3‏ با هم همپوشانی 
شدند و نقشه شدت فرسایش منطقه به دست آمد؛ بنابراین و بر اساس رابطه ۱۳ و با همپوشانی لایه‌های 
اطلاعاتی. نقشه شدت فرسایش حوضه (شکل ag (VA‏ گردید. لازم به ذکر است ضریب شدت فرسایش 
(Z)‏ برای هر یک از واحدهای همگن به‌صورت جداگانه محاسبه گردید (جدول A‏ و میانگین OF‏ برای کل 


حوضه در نظر گرفته شد. بر اساس اطلاعات به دست آمده از جدول A‏ ضریب استفاده از زمین حداقل ۰/۲۰ 


۱۶۰۳ تابستان‎ Y شماره‎ Y نشریه جغرافیا و مخاطرات محیطیء جلد‎ YF. 


و حداکثر ۰/۰ می‌باشد و متوسط ol‏ در مدل EPM‏ ۰/۶۰ برآورد شده است. متوسط ضریب حساسیت 
سنگ و خاک به فرسایش CY)‏ ۰/۷۱ می‌باشد که کمترین آن ۰/۵۰ و بیشترین آن ۱/۲۰ بر آورد شده است. 
ضریب فرسایش مشاهده‌ای یا حوضه منطقه که حداقل آن ۰/۲۰ و حداکثر Ol‏ ۰/۷۰ برآورد شده است. 
ضریب فرسایش Z)‏ در این منطقه به‌طور متوسط در کل حوضه برابر ۰/۸۹ می‌باشد که بر اساس جدول ۳ 
در محدوده فرسایش 11 یا شدید قرار دارد. جدول ۸ ضرایب فرسایش حوضه و ضریب شدت فرسایش را 
در واحدهای همگن به تفکیک مساحت نشان داده است. نتایج حاصل از بررسی نقشه شدت فرسایش OUS‏ 
می‌دهد که بیشترین فرسایش در قسمت غرب و شرق حوضه قرار دارد. چون این منطقه از سازندهای ناپایدار 
و مستعد فرسایش تشکیل شده است. شیب obj‏ منطقه به‌اضافه بارش باران و تغییر کاربری اراضی در این 
محدود بسیاری از عرصه‌های ol‏ را با بحران فرسایش خاک مواجه کرده است. مقدار Z‏ برای این مناطق بیش 
از Y‏ است که نشان‌دهنده فرسایش‌پذیری شدید این ناحیه می‌باشد(شکل (VA‏ 


حدول ^- داده‌های Ce e‏ محاسبه cA‏ فرسایش خاک در واحدهای همگن با EPM Ju.‏ 


Table 8- Necessary data to calculate soil erosion intensity in homogeneous units with EPM 


model 

ale‏ ضریب ضریب ضریب ‏ شیب 0 ضریب کلاس مساحت 

واحد استفاده Cagle‏ عامل متوسط شیب شدت فرسایش (هکتار) 

همگن اززمین سنگ و خاک فرسایش (در صد) فرسایش 

ha(area) Erosion Z 1 slope Y y Xa PN 

intensity 

MYY duas UNE Ye 3 T ۰/۲ UA \ 
۹0۰/4 ETE Uv MAU ۷ TZ Fat ۰۳۵ ۲ 
۱۷ متوسط‎ Y ۳/۰۰ 4 ۰/0۰ VAL AME Y 
AY 0/AY شد ید‎ 2 Y/M ۷ ۰/۳۰ vas ۰/0۰ 3 
۱۶۵ خیلی شدید‎ ۱/۸۰ ۳/۹ E AL VY ۰/۰ ۵ 
۸3 متوسط‎ +/00 ۳/۳۹۶ ۷ ۰/10۵ ۰۹۰ NS a 
۱۸ متوسط‎ MAT AU ۳ EE Y NS ۷ 
(۱۳۸۹ شد ید‎ ۰/۸۱ 1۹۹ ۳۲ AE NV ۰/۳۲ ۸ 
2۱۰۹/۳ خیلی شدید‎ ۱/۷ 0/0 YY YY ۰/۳۲ UNS 4 
YYVV/AY متوسط‎ ۰/۵ ۱/۰ ۱ AL ۰/۳۲ ۳۵ E 
16/۰۷ متو سط‎ V/A \/£\ Y ۰/0۰ ۱/۳۲ Ns \\ 
ON NA خیلی شدید‎ EA \/2\ Y */00 MY ۰/۹۰ yy 


4/۲7 متو سط‎ Y ۱/۱ Y JN 1/08 A ۱۳ 
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۱۱۹۱/۹۱ شدید‎ ۰۹۹ 11۷ Y. ۰/0۰ 1/07 EE M 
£*/0^ خیلی شدید‎ ۱/۰۰ Y/YY ۵ +/00 ۱/۲۰ AN ۱۵ 
58370 ۰/۸۹ ۲/۸۰ ۵۳ ۰/۸ AY ۰/۳۷ میانگین‎ 


۵ و ۱۷ شان osla‏ شده است. همچنین در شکل ۱۸ نقشه شدت فرسایش در حوضه چم‌گردلان Obss‏ 
oala‏ شده است که بر اساس GT‏ کمترین مقدار فرسایش در قسمت‌های مرکزی و هموارتر منطقه است و بیشتر 


زیر کشت آبی و دیم و باغ‌کاری در مناطق با شیب کم قرار می‌گیرد که دارای فرسایش‌پذیری کمتری می- 
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شکل ۱۶- ضریب فرسایش حوضه(۳) شکل ۱۵- ضریب حساسیت سنگ و خاک به فرسایش(۷) 
Fig.14. basin erosion coefficient (¥)‏ 


Fig.15. Coefficient of sensitivity of rock and soil to erosion (Y) 


نشریه جغرافیا و مخاطرات محیطی» جلد ۰۱۳ شماره ۲ تابستان ۱۶۰۳ 
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(Xa) ie ضریب استفاده از‎ -N شکل‎ 
Fig.16. Land use coefficient (Xa) 
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۲12.17: Average slope coefficient of the basin (I) 
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شکل ۱۸- ضریب شدت فرسایش(7) 
Fig.18. Erosion intensity coefficient (Z)‏ 


گردیده و با همپوشانی نقشه سالانه فرسایش ویژه و فرسایش کل در حوضه به دست آمد(شکل‌های ۱٩‏ و 
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Fig.19. Special erosion 


شکل ۲۰- فرسایش کل 


Fig.20. total erosion 


درنهایت از تلفیق نتایج به دست آمده» فرسایش ویژه فرسایش JS‏ رسوب ویژه و رسوب کل در واحدهای 
همگن به دست آمد (جدول A‏ بر این اساس» میزان فرسایش ویژه و فرسایش کل به ترتیب معادل ۱۵۰۲ 
مترمکعب در کیلومترمربع در سال و ۸۷۷۱۷۶ مترمکعب در سال است. همچنین ميزان رسوب ویژه و رسوب 
کل در حوضه آبخیز چم‌گردلان به ترتیب ۲۹۵۹ مترمکعب در کیلومترمربع در سال و ۱۱۹۱۰۰۲ مترمکعب در 
سال است. همچنین در شکل‌های ۲۱ ۲۲ و ۲۳ به ترتیب نقشه‌های ضریب رسوب. رسوب ویژه و رسوب 
کل نشان داده ott‏ است. با توجه به جدول شماره ٩‏ میزان فرسایش ویژه بسته به ميزان ضریب شدت 
فرسایش در واحدهای همگن تغییر می‌کند. این تغییرپذیری» نتیجه کاربری‌های مختلف حوضه و همچنین 
مااع مات و تست oldies‏ 

مقدار رسوب ویژه برای واحدهای همگن و کل حوضه بر حسب مترمکعب در کیلومترمربع در سال به دست 
آمد. سپس با در نظر گرفتن وزن مخصوص ظاهری خاک VEN)‏ مقدار رسوب ویژه بر حسب تن در هکتار 
در سال محاسبه شد و مقدار YVES Of‏ تن در هکتار در سال بر آورد گردید (جدول 4). این در Se‏ است 
که مقدار رسوب اندازه‌گیری شده در ایستگاه‌های هیدرومتری و رسوب سنجی برابر ۱۳/۵۸ تن در هکتار در 
سال می‌باشد؛ به‌عبارت‌دیگر مقایسه نتایج ایستگاه هیدرومتری و رسوب سنجی در خروجی حوضه با مدل 
EPM‏ نشان‌دهنده هم‌خوانی VY‏ درصدی رسوبی خروجی این مدل است. در جدول ۱۰ مقادیر مشاهده شده 


و برآورد شده در حوضه آبخیز چم گردان استان ایلام را مشاهده می‌فرمایید. 


۱۶۰۳ شماره ۲. تابستان‎ VY جغرافیا و مخاطرات محیطی. جلد‎ 4 pia 


۳۶۴ 


حدول 4— فرسایش و رسوب محاسبه شده برای واحدهای همگن و کل حوضه با مدل EPM‏ 


Table 9- Erosion and sedimentation calculated for homogeneous units and the entire basin 


with the EPM model 
5. ^ NC شماره واحد فرسایش ویژه مساحت فرسایش کل ضریب رسوبدهی‎ 
m/year nr /year/ 2 
( y ) ( y ) مترمکعب در سال) در واحدهای همگن‎ ( (km?) همگن (مترمکعب‎ 
در کیلومترمربع در سال)‎ 
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شکل ۲۱- ضریب رسوبدهی(181) 


Fig.21. Sedimentation coefficient (RU) 


شکل ۲۲- رسوب ویژه 
Fig.22. Special sediment‏ 
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۲12.23. Total sediment 


در جدول ۱۰ polis‏ مشاهده شده و برآورد شده در حوضه آبخیز چم‌گردان استان ایلام را مشاهده 
می‌فر مایید. نتایج به دست achat‏ نشان می دهد که ميزان رسوب برآوردی مدل RUSLE‏ ۱/۳۷ تن در هکتار 
در سال و در روش EPM‏ ۹ تن در هکتار در سال برآورد 5 19 که با توجه به نزدیک‌تر بودن ET‏ 
برآوردی RUSLE‏ به مقدار ow)‏ — مشاهده‌ای ۱۳۹/۵۸ تن در هکتار در سال» کارایی مدل RUSLE‏ در 
حوضه از روش RUSLE‏ بهره گرفت. 

حدول ۰- مقایسه polis‏ مشاهده شده و برآورد شده در حوضه آبخیز چمگردان استان ایلام 


Table 10- Comparison of the observed and estimated values in Chamgardan aquifer basin 
of Ilam province 


مقدار مشاهده شده مقدار برآورد شده با مدل RUSLE‏ مقدار برآورد شده با مدل EPM‏ 
(تن در هکتار درسال) (تن در هکتار درسال) (تن در هکتار درسال) 
YAS VV Y'V6A^‏ 
نتیجه گیری 


اندازه گیری مقدار فرسایش و رسوب به دلایل فنی. حفاظتی و اقتصادی در تمامی نقاط امکان‌پذیر نیست. 
ازاین‌رو مناسب‌ترین راهکار برآورد مقدار فرسایش و تولید رسوب با مدل‌های تجربی است که لازمه آن 
آگاهی از مکانیسم‌های فرسایش و عوامل مؤثر بر آن است. از طرفی انتخاب مدل مناسب برای هر منطقه 
مستلزم ارزیابی دقت مدل‌های مختلف از طریق مقایسه نتایج حاصل از به‌کارگیری آن‌ها به ضریب رسوبدهی 


واقعی در یک حوضه آبخیز می‌باشد. بر این اساس» در این تحقیق, میزان فرسایش و رسوب حوضه آبخیز 
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چم‌گردلان ایلام با استفاده از دو مدل EPM‏ و RUSLE‏ موردبررسی قرار گرفت. از طرفی میزان رسوب 
مشاهده‌ای با استفاده از روش منحنی سنجه رسوب (حدواسط دسته‌ها) برآورد گردید و درنهایت موثرترین و 
کاراترین روش معرفی گردید. 

در این مطالعه به‌منظور بررسی میزان فرسایش و برآورد مقدار رسوب در حوضه آبخیز چم‌گردلان استان 
ایلام» پس از امتیازدهی و تلفیق عوامل ij‏ به‌صورت لایه‌های رقومی در محیط نرم‌افزار ArcGIS‏ با 
استفاده از مدل‌های تجربی نقشه میزان فرسایش و رسوب حوضه موردنظر تعیین گردید. 

نتایج برآورد معیارهای موردبررسی در روش RUSLE‏ در حوضه آبخیز موردمطالعه نشان داد. فاکتور 
فرسایش پذیری از ۱۰/۲۸ ۰/۱۱ و معیار پوشش گیاهی از ۰/۳۵ تا ۰/۷۰ متغیر بوده که نشان‌دهنده تفاوت در 
ویژگی‌های خاک و پوشش گیاهی در مناطق مختلف حوضه است. نقشه فاکتور پوشش گیاهی نشان می‌دهد 
این پارامتر در مناطق بدون پوشش گیاهی حداکثر است و کمترین مقدار ol‏ در مناطقی است که پوشش 
متراکم وجود دارد. نتایج حاصل از فاکتور فرسایندگی در حوضه چم‌گردلان نشان داد. میزان این فاکتور از 
۲ در پایین‌دست تا YA‏ مگاژول در میلی‌متر در سال در هکتار در بالادست حوضه در نوسان است و میزان 
فاکتور فرسایندگی باران در نقطه حروجی حوضه به حداقل می‌رسد. Miguel, Samuel-Rosa, Simao‏ 
Bahrawi et al.(2016) ‘Dennis Dalmolin, Arajo Pedron, & Moura Bueno (2011)‏ و Rezaee,‏ 
Faridi, Ghorbani, & Kazemi (2018)‏ معتقدند افزایش ارتفاع از سطح درا سبب افزایش فاکتور 
فرسایندگی باران می‌گردد. چون که میزان بارش زیاد می‌شود. فرسایندگی باران در شمال‌شرقی و جنوب‌غربی 
حوضه موردمطالعه که ارتفاع بالاتری داشتند بیش‌تر بوده که با رنگ قرمز نشان داده شده است یا 
به‌عبارت‌دیگر کاهش ارتفاع و باران رابطه مستقیمی با کاهش فاکتور R‏ دارد. بررسی همبستگی آماری بین 
فاکتورهای RUSLE‏ و مقدار فرسایش سالیانه خاک در منطقه موردمطالعه نیز نشان داد که فاکتور LS‏ با 
ضریب همبستگی ۰/۸۱ بیشترین اهمیت را در برآورد فرسایش سالانه خاک در مدل RUSLE‏ دارد. 

در انجام محاسبات مربوط به مدل EPM‏ میانگین ضرایب فرسایش حوضه (P)‏ استفاده از زمین Xa)‏ 
حساسیت سنگ و خاک به فرسایش (Y)‏ و شیب متوسط (D)‏ محاسبه و با مقایسه با جدول استاندارده ميزان 
فرسایش و کلاس فرسایش برای هرکدام از واحدهای همگن و درنهایت در کل حوضه برآورد گردید. 
جنانچه در جدول مربوطه مشاهده می‌فرمایید. ضریب شدت فرسایش E)‏ به‌جز ٩‏ واحد از ۱۵ واحد در 
کلاس شدت فرسایش شدید تا خیلی شدید قرارداشته و میانگین Ol‏ در کل حوضه ۰/۸۹ برآورد گردید. 
همچنین میزان فرسایش ویژه و رسوب ویژه کل حوضه به ترتیب ۱۵۰۲/۷۷ و ۲۹۵۳/۷۸ مترمکعب در 
کیلومترمربع در سال برآورد گردید. نتایج حاصل از کاربرد EPM‏ در برآورد فرسایش OLS‏ داد که میزان VO‏ 


در صد از مساحت حوضه در کلاس فرسایشی شدید (زیاد) و خیلی شدید GU»)‏ دارد (شرقی و میانی 
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حوضه) و حدود ۳۵ درصد حوضه در کلاس فرسایش کم تا متوسط است. با استفاده از روش فوق. مقدار 
رسوب سالانه YVES‏ تن در هکتار در سال برآورد گردید. 

بر اساس نتایج fob‏ از روش RUSLE‏ بیشترین وسعت حوضه مربوط به کلاس فرسایشی خیلی کم کم و 
متوسط است که عموماً در کل حوضه پراکنش دارد و کمترین مساحت حوضه در کلاس فرسایشی زياد تا 
he‏ زیاد (۲۰ درصد) قرار می گیرد و بر این اساس مقدار رسوب سالانه در این روش ۱3/۳۷ تن در JUS‏ 
در سال برآورد گردید. بنابراین» نتایج حاصل از EPM‏ غالب فرسایش را در کلاس شدید تا خیلی شدید و 
مدل RUSLE‏ در کلاس خیلی‌کم تا متوسط برآورد نموده‌اند» بدین معنی که مدل EPM‏ دچار پیش‌برآورد 
و مدل RUSLE‏ دچار کم‌برآورده شده است. نتایج حاصل از بررسی نقشه فرسایش در هر دو مدل OLS‏ 
eae caca‏ جنس Bb cx‏ عرش ترا هرگ عون ای مان Jie cia‏ 
و مستعد فرسایش تشکیل شده است. شیب زیاد منطقه به اضافه بارش OL‏ و تغییر کاربری اراضی در این 
محدود بسیاری از عرصه‌های آن را با بحران فرسایش خاک مواجه کرده است. مقدار Z‏ برای این مناطق پیش 
از ۱ است که نشان‌دهنده فرسایش‌پذیری شدید این منطقه می‌باشد. 

یکی از عوامل موثر بر فرسایش و تولید رسوب. رژیم بارندگی و اقلیم حوضه است که با بررسی آن مشخص 
شد که قسمت عمده ریزش‌های جوی در حوضه آبخیز چم‌گردلان از نوع باران است و در فصل زمستان با 
برف همراه می‌شود ولیکن بارندگی به‌عصوص با شدت‌های بالا و زمان کې باعث جریان سریع و شدت 
رواناب شده که با خود مقدار زیادی از خاک سطحی را تخریب و به‌صورت رسوب و بار معلق به مناطق 
پایین‌دست و خروجی حوضه حمل می‌کند. 

از دیگر عوامل موثر بر فرسایش و تولید رسوب و ایجاد سیلاب حوضه. خحصوصیات فیزیوگرافی و 
توپوگرافی ازجمله شیب و ارتفاع است. شیب حوضه آبخیز چم گردلان در اکثر قسمت‌های آن چشمگیر 
است به‌طوری که در بیشتر قسمت‌های این حوضه شیب بیشتر از ۲۰ درصد است که در سرعت بخشیدن به 
فرسایش و پارامترهای مؤثر در آن ازجمله رواناب NS SE‏ است. به‌طورکلی در حوضه‌های با شیب ملایم به 
دلیل اینکه قدرت آب جاری کاهش پیدا می‌کند. بنابراین اشکال مختلف فرسایش نیز کمتر دیده می‌شود؛ 
db‏ که در حوضه آبخیز چم‌گردلان به دلیل کوهستانی بودن قسمت‌های زیادی از آن. عکس این حالت 
اتفاق می‌افتد. 

پوشش گیاهی و کاربری اراضی یکی دیگر از عوامل مؤثر در فرسایش و رسوب حوضه به شمار می‌رود. 
طبیعی است مناطقی که دارای پوشش گیاهی جنگلی و مرتعی مناسب باشد. کمتر در معرض تخریب و 
فرسایش قرار می گیرند. در حوضه آبخیز چم گردلان که حدود ۷۰ درصد از مساحت حوضه را جنگل و مرتع 
با تراکم‌های مختلف پوشانده است. انتظار می‌رود. وضعیت کنترل فرسایش در سطح حوضه آبخیز مناسب و 
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میزان ضریب فرسایش در حد کم باشد؛ لیکن به دلایل مختلف ازجمله چرای مفرط و دائم مرتع توسط دام- 
های منطقه عدم رعایت تعادل دام و مرتع توسط دامداران. اعمال روش‌های غلط شخم اراضیء قطع بی‌رویه 
درختان» درختچه‌ها و بوته‌هاء عدم رعایت اصول بهره‌برداری صحیح از اراضی. تبدیل جنگل‌ها و مراتع به 
اراضی کشاورزی» SE‏ و زمین‌شناسی و SE‏ عوامل اقلیمی Ge pete‏ پراکنش و شدت بارندگی» باعث 
شده است که وضعیت فرسایش حوضه آبخیز در حد شدید و بسیار شدید قرار pF‏ 2 لذا اعمال روش‌های 


مختلف کنترل فرسایش در سطح حوضه ضروری است. 

نتایج به دست آمده OUS‏ می‌دهد میزان رسوب برآوردی در روش RUSLE‏ ۱۱/۳۷ تن در هکتار در سال و 
در روش ۲۱/4٩ EPM‏ تن در هکتار در سال برآورد گردید که با توجه به نزدیک‌تر بودن مقدار برآوردی 
RUSLE‏ 4 مقدار مشاهده‌ای ۱7/۵۸ تن در هکتار در سال» کارایی مدل RUSLE‏ در برآورد رسوب منطقه 
موردمطالعه بیشتر است. محققینی چون Vaezi, Bahrami, ‘Miguel et al., (2011)*Bahrawi et al.(2016)‏ 
Teimourri & (2018) Rezaee, Faridi, Ghorbani, & Kazemi, ‘Sadeghi, & Mahdian (2010)‏ 
)2018( 2221015020 کارایی بالای مدل RUSLE‏ را ast‏ می‌نمایند.(1995) Abedini, ¢ Rafahi & Nemeti‏ 
Khodabakhsh, Mohammadi, Rafie, & Bozorg Zadeh, (2010) : Shabrang, & Esmaeli (2014).‏ و 
Sadog, Hoseinzadeh, & Azadi (2015)‏ کارایی نامناسب مدل EPM‏ را در برآورد فرسایش و رسوب در 


مطالعات تفصیلی ASL‏ نموده‌اند. 


References 


Abedini, M., Shabrang, S., & Esmaeli, A. (2014). Review on Soil Erosion and Sedimentation in 
Meshkin Chai Catchment Area by EPM Method. Geography and Development, 11(30), 
87-100. [In Persian] https://doi.org/10.2211 1/gdij.2014.245 

Amiri, F., Arzani, H., Farahpour, M., Chaeichi, M.R., & Khajeh Aldin, S.J.A.D. (2009). 
Efficiency of Mpsiac and EPM Models for Assessment Soil Erosion in Range 
Suitability. Rangeland, 3(1), 138-154. [In Persian] https://sid.ir/paper/136360/en 

Ansari Lari, A., & Ansari, M. (2017). Evaluation of soil erosion risk and sedimentation 
potential by using EPM model in Gabric basin- SE Hormozgan- Iran. Journal of 
Natural Environmental Hazards, 6(11), 1-14. [In Persian] 
https://doi.org/10.22111/jneh.2017.2965 

Asgari, S.A., Servati, M.R., & Jafari, M.R. (2008). Estimating the Oil Erosion and Production 
of Sediment in the Basin of the Ilam Dam by using MPSIAC Model. Physical 
Geography Research, 40(64), 29-35. [In Persian] https://www.sid.ir/paper/5427/en 

Azami, A. (2001). Comparison of Various Methods for Estimating Sediment Yield in Ilam Dam 
Watershed. (published M.sc. thesis). Gorgan University. (in persian). 

Azimi Sardari, M.R., Bazrafshan, O.A., Panagupulus, T., & Rafiei Sardoei, E. (2019). Current 
and Future Assessment of Soil Erosion in the Catchment Area of Esteghlal Dam in 


سال سیزدهم آرخی و همکاران. مقایسه مدل‌های RUSLE‏ و EPM‏ در برآورد فرسایش و تولید رسوب با استفاده از ... ۳۶۹ 


Minab Using RUSLE-3D Model and Climate Change Scenarios. Desert Management, 7 
(14), 132-119. [In Persian]. 
https://doi.org/10.22034/jdmal.2020.38480. 

Bahrawi, A.J., Elhag, M., Aldhebiani, A.Y., Galal, H.K., Hegazy, A.K., & Alghailani, E. (2016). 
Soil Erosion Estimation Using Remote Sensing Techniques in Wadi Yalamlam Basin, 
Saudi Arabia. | Advances in Materials Science and Engineering. 
https://doi.org/10.1155/2016/9585962. 

Bartsch, K.P., Van Miegroet, H., Boettinger, J., & Dobrwolski, J.P. (2002). Using Empirical 
Erosion Models and GIS to Determine Erosion Risk at Camp Williams. Journal of Soil 
and Water Conservation, 57(1), 29-37. https://www.jswconline.org/content/5 7/1/29/tab- 
references. 

Benkobi, L., Trlica, M.J., & Smith, J.L. (1994). Evaluation of a Refined Surface Cover 
Subfactor for Use in RUSLE. Journal of Range Manage 47(1),74—78. 

https://doi.org/10.2307/4002845 

Biesemans, J., Meirvenne, M.V., & Gabriels, D. (2000). Extending the RUSLE with the Monte 
Carlo Error Propagation Technique to Predict Long-Term Average Off-Site Sediment 
Accumulation. Journal of Soil and Water Conservation, 55(1), 35-42. 
https://www.jswconline.org/content/55/1/35 

Bagio, á., Bertol, I., Wolschick, N.H., Schneiders, D., & Aparecida, M. (2017). Water Erosion 
in Different Slope Lengths on Bare Soil. Revista Brasileira de Ciéncia do Solo. 

DOI: 10.1590/18069657rbcs20160132 

Boggs, G., Devonport, C., Evans, K., & Puig, P. (2001). GIS-Based Rapid Assessment of 
Erosion Risk in a Small Catchment in the Wet/Dry Tropics of Australia. Land 
Degradation and Development, 12(5), 417—434. https://doi.org/10.1002/ldr.457 

Boyce, R. C. (1975). Sediment Routing with Sediment Delivery Ratios. In: Present and 
Prospective Technology for Predicting Sediment Yields and Sources, US Dept. Agric. 
Publ. ARS-S-40, 61-65. 

Deore, S.J. (2006). Prioritization of Micro-watersheds of Upper Bhama Basin on 
the Basis of Soil Erosion Risk Using Remote Sensing and GIS Technology. 
(published doctoral. dissertation). University of Pune, Department of Geography. 

Elirehema, Y. (2001). Soil Water Erosion Modeling in Selected Watersheds in Southern Spain. 
IFA, ITC, Enschede. https://doi.org/10.1002/1dr.457 

Ferro, V., Giordano, G., & Lovino, M. (1991). Isoerosivity and Erosion Risk Map for Sicily. 
Hydrological Science Journal, 36(6),549—564. 
https://doi.org/10.1080/02626669 109492543. 

Gavrilovic, Z. (1988, May). The Use of an Empirical Method (Erosion Potential Method) for 
Calculating Sediment Production and Transportation in Unstudied or Torrential streams. 
Paper presented at the International Conference of River Regime, Wallingford. 
https://www.scirp.org/reference/referencespapers?referenceid=1840329. 

Haan, C.T., Barfield, B.J., & Hayes, J.C. (1994). Design Hydrology and Sedimentology for 
Small Catchments. San Diego, Academic Press. 
https://www.scirp.org/reference/referencespapers ?referenceid=1687737 

Khodabakhsh, S., Mohammadi, A., Rafie, B., & Bozorg Zadeh, I. (2010). Comparison of 
erosion and sediment yield estimation in Sezar sub-basin (Dez drainage basin) by 
MPSIAC and EPM empirical methods, using GIS. Iranian Journal of Geology, 3(12), 
51-61. [In Persian] https://www.magiran.com/p886221 

Lal, R. (1990). Soil Erosion in the Tropics: Principles and Management. New York, McGraw- 
Hill. https://www.scirp.org/reference/referencespapers?referenceid=2 193850 


۳۷۰ نشریه جغرافیا و مخاطرات محیطیء جلد Y‏ شماره Y‏ تابستان ۱۶۰۳ 


Lin, C.Y. (1997). A Study on the Width and Placement of Vegetated Buffer Strips in a 
Mudstone-Distributed Watershed. Journal of China. Soil Water Conserve, 29 (3), 250- 
266 [In Chinese with English abstract]. 

McCool, D.K., Brown, L.C., & Foster, G.R. (1987). Revised Slope Steepness Factor for the 
Universal Soil Loss Equation. Transactions of American Society of Agricultural 
Engineers, 30, 1387-1396. https://doi.org/10.13031/2013.30576. 

Miguel, P.A., Samuel-Rosa, R., Simao Dennis Dalmolin, F., Arajo Pedron, J., & Moura Bueno, 
A. (2011). The USLE Model for Estimating Soil Erosion in Complex Topography 
Areas. Annals XV Brazilian Symposium on Remote Sensing, (SBSR), Brasil, 85, 9227- 
9230. 

Milevski, I. (2008, July). Estimation of Soil Erosion Risk in the Upper Part of Bregalnica 
Watershed-Republic of Macedonia, Based on Digital Elevation Model and Satellite 
Imagery. Paper presented at the 5th International Conference on Geographic 
Information Systems (ICGIS-)(Turkey), Istanbul,( pp.351-358). 
https://www.researchgate.net/publication/265012436 

Mokhtari, L.G., Shafiei, N,. & Rahmani, A. (2019). Estimation of Soil Erosion Using RUSLE 

Model, A Case Study of Noorabad Mamasani Basin. Hydrogeomorphology, 5(17), 1-21. 

[In Persian). https://hyd.tabrizu.ac.ir/article 8601. en.html?lang-en 

Moore, I.D., & Wilson, J.P. (1993). Length-Slope Factors for the Revised Universal Soil Loss 

Equation: Simplified Method of Estimation. Journal of Soil and Water Conservation 

47(4),423-428. https://www.scirp.org/reference/referencespapers?referenceid=674926. 

Moradi, M., Ghonchepour, D., Nohegar, A., & Mahmoudi Nejad, V. (2012). A Comparison of 
the MPSIAC and EPM Models for Estimating Erosion and Sediment in the Poorahmadi 
Catchment. Environmental Erosion Researches, 1(4), 54-68. [In Persian] 
https://www.sid.ir/paper/223084/en 

Morgan, R.P.C. (1995). Soil Erosion and Conservation. London, Longman. 
https://www.scirp.org/reference/referencespapers?referenceid=117218 

Noori, H., Siadatmousavi, S. M., & Mojaradi, B. (2017). Estimation of Soil Erosion and 
Sediment Yield Based on Satellite Imagery Using the GIS Technique and Employing 
Erosion Potential Method (EPM) on Watersheds (Case study: the Dez Watershed). 
Water Resources Engineering, 10(32), 45-54. 
https://dorl.net/dor/20.1001.1.20086377.1396.10.32.6.2 

Nouwakpo, S.K., Williams, C.J., Al-Hamdan, O.Z., Weltz, M.A., Pierson, F. & Nearing, M. 
(2016). A Review of Concentrated Flow Erosion Processes on Rangelands: 
Fundamental Understanding and Knowledge Gaps. /nternational Soil and Water 
Conservation Research, 4(2), 75-86. https://doi.org/10.1016/j.iswer.2016.05.003 

Rafahi, H., & Nemeti, M. (1995). Application of EPM Method in Studying Rate of Erodibility 
and Sedimentation in Alamut-Rud Watershed. /ranian Journal of Agricultur Science, 
26(1), 33-46. [In Persian] 
https://journal.ut.ac.ir/article 16105 60140101923 5 1 

Renard, K.G., & Freimund, J.R. (1994). Using Monthly Precipitation Data to Estimate the R 
Factor in the Revised USLE. Journal of Hydrology, 157(1-4), 287—306. 
https://doi.org/10.1016/0022-1694(94)90110-4. 

Renard, K.G., Foster, G.R., Weesies, G.A., McCool, D.K., & Yoder, D.C. (1997). Predicting 
Soil Erosion by Water: A Guide to Conservation Planning with the Revised Universal 
Soil Loss Equation (RUSLE). Agriculture Handbook, 703, USDA-ARS. 
https://www3.epa.gov/npdes/pubs/ruslech2.pdf 

Rezaee, P., Faridi, P., Ghorbani, M., & Kazemi, M. (2018). Estimation of Soil Erosion Using 
the RUSLE Model and Identification of its Most Effective Factor in the Gabrik 


سال سیزدهم آرخی و همکاران. مقایسه مدل‌های RUSLE‏ و EPM‏ در برآورد فرسایش و تولید رسوب با استفاده از ... ۳۷/۱ 


Watershed-South-Eastern Hormozgan Province Quantitative Geomorphological 
Research, 3(1), 97-113. [In Persian] 
https://dorl.net/dor/20.1001.1.22519424.1393.3.1.7.3 

Sadog, H., Hoseinzadeh, M.M., & Azadi, F. (2015). Determining the Erosion in Kahman 
Drainage Basin Using EPM, BLM and Fargas Models. Hydrogeomorphology, 2(2), 
137-154. [In Persian] https://dorl.net/dor/20.1001.1.23833254.1394.2.2.8.7 

Shahbazi, K. (1999). Estimation of Erosion and Sediment by Qualitative Geomorphology 
(Homogenous Units) and EPM Methods and Comparing them with Outlet Sediment 
Data in Ilam Dam Basin. (published MS.c thesis). Tarbiat Modarres University, School 
of Natural Resources. [In Persian]. 

Tazioli, A. (2009). Evaluation of Erosion in Equipped Basins: Preliminary Results of a 
Comparison Between the Gavrilovic Model and Direct Measurements of Sediment 
Transport. Environmental Geology, 56(5), 825-831. 
https://doi.org/10.1007/s00254-007-1183-y 

Teimourri, F., & Bazrafshan, O. (2018). Comparison of the Efficiency the RUSLE and EPM 
Models in Estimating Soil Erosion and Sediment Yield According to the Sediment 
Rating curve (case study :Jask-Gabrik watershed). New Findings in Applied 
Geology,12(23), 50-62. [In Persian] 
https://journals.basu.ac.ir/article 2278 101b5eaf999bfabc46cbcf8240155cle.pdf 

USDA. (1975). Sediment Sources, Yields, and Delivery Ratios. National Engineering 
Handbook, Section 3 Sedimentation. 
https://irrigationtoolbox.com/NEH/Part632 SedimentationGeology/NEH3 06.PDF 

Vaezi, A.R., Abasi, M., & Abd Haji Maleki, K.H. (2017). Evaluation of RUSLE Model 
Combined with Remote Sensing and GIS in Small Drainage Areas in the Semi-Arid 
Region, Northwestern Iran. /ranian Journal of Watershed Management Science and 
Engineering, 11 (38), 1-10. [In Persian]. 
https://www.sid.ir/paper/134808/en#downloadbottom. 

Vaezi, A.R., Bahrami, H.A., Sadeghi, S.H.R., & Mahdian, M.H. (2010). Spatial Variability of 
Soil Erodibility Factor (K) of the USLE in North West of Iran. Journal of Agricultural 
Science and Technology, 12, 241-252. [In Persian]. 
http://dorl.net/dor/20.1001.1.16807073.2010.12.2.11.9 

Vanoni, V.A. (1975). Sedimentation Engineering, Manual and Report No. 54. American Society 
of Civil Engineers, New York, N.Y.  https://webstore.ansi.org/preview- 
pages/ASCE/preview 9780784471340.pdf. 

Wischmeier, W.H., & Smith, D.D. (1978). Predicting Rainfall Erosion. Losses: A Guide to 
Conservation Planning. Agriculture Handbook, 537. US Department of Agriculture, 
Washington, DC,. 
https://www.ars.usda.gov/ARSUserFiles/60600505/RUSLE/AH_537%20Predicting%2 
ORainfall%20So0i1%20Losses.pdf 

Yu, B., & Rosewell, C.J. (1996). A Robust Estimator of the R Factor for the Universal Soil Loss 
Equation. American Society of Agricultural Engineers, 39(2), 559-561. 
https://typeset.io/pdf/a-robust-estimators-of-the-r-reaction-for-the-universal-soil- 
7rxy 1i22im.pdf 

Yue-Qing, X., Shao, X-M., & Kong, X-B. (2008). Adapting the RUSLE and GIS to Model Soil 
Erosion Risk in a Mountain Karst Watershed, Guizhou Province, China, Environment 
Monitoring and Assessment, 141, 275-286. https://doi.org/10.1007/s10661-007-9894-9. 


